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Abstract:  A  study  of  blade-surface  asperities  with  respect  to  their  effect  on 
cavitation  onset  on  marine  propellers  is  carried  out  to  investigate 
possible  quality  control  criteria  for  propeller  models  intended  for 
cavitation  inception  tests.  The  0B0  propeller  is  used  as  an 
example  to  test  possible  calculation  procedures.  ~u  asperities 
consisted  of  two-dimensional  isolated  triangular  roughness  elements, 
distributed  roughness  and  larger  scale  two-dimensional  blauc-surface 
bumps  or  waves.  It  was  found  that  useful  quality  control  criteria 
cannot  be  formulated  for  triangular  roughness  elements  because  the 
basic  data  which  would  permit  their  development  are  not  available. 
However,  it  does  appear  feasible  to  develop  quality  control  criteria 
for  the  other  two  types  of  blade  asperities  and  computation  procedures 
which  can  lead  to  the  formulation  of  blade-surface  fabrication 
specifications  are  given  for  them. 
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NOMENCLATURE 

profile  chord  length 

1  2 

pressure  coefficient,  (p-  p  )/  —  p  V 

o  Z.  o 

minimum  pressure  coefficient 
pressure  coefficient  on  blade  surface 

section  lift  coefficient 

section  lift  coefficient  at  ideal  attack  angle  (design  lift 
coefficient) 

relative  error  in  cavitation  number  caused  by  blade  surface 
asperity,  (acR  -  ac£.)/acs 

maximum  camber  of  blade  profile 

maximum  camber  of  profile  or  maximum  height  of  surface  asperity 
camber  correction  factor 

Reynolds  number  based  c  blade  section  chord 
free  stream  velocity 

skin  friction  velocity,  /i^/p 

distance  along  chord  line 
distance  normal  to  chord  line 

angle  of  attack  measured  with  respect  to  the  chord  line 
ideal  or  "design"  angle  of  attack 

mean  relative  inflow  angle  measured  with  respect  to  propeller 
plane 

variations  in  propeller  inflow  angle  caused  by  hull-induced 
inflow  distortions 


1980 


boundary  layer  thickness 
displacement  thickness 


-5- 


August  20,  1980 
F’^tpjk 


momentum  thickness 


viscous  sublayer  thickness 


dimensionless  bump  height  h/£ 
ordinate  of  bump 
kinematic  viscosity 
fluid  density 
length  of  limb 

1  2 

cavitation  number,  (p  -  p  )/  t  p  V_ 

o  v  /.  c 

inception  cavitation  number  on  surface  of  a  smooth  blade 


incipient  cavitation  number  on  roughened  body  or  on  body  with 
asperity 

cavitation  number  of  an  asperity  on  a  flat  plate 
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INTRODUCTION 

Experience  vith  model  propeller  tests  in  cavitation  tunnels  and 
towing  tanks  seems  to  suggest  that  quality  control  requirements  with 
respect  to  blade-surface  contours  and  roughness  may  be  more  severe  for 
cavitation  onset  measurements  than  is  the  case  for  overall  propulsive 
performance  evaluation.  Therefore,  the  16th  ITTC  Cavitation  Committee 
has  been  assigned  the  task  of  determining  whether  or  not  sufficient  basic 
data  are  now  available  to  enable  the  formulation  of  quality  control 
recommendations  which  member  organizations  can  use  when  fabricating  propeller 
models  for  cavitation  testing.  This  report  documents  one  aspect  of  a  more 
comprehensive  review  of  the  subject  carried  out  by  the  committee  upon  the 
effects  of  blade  surface  roughness  and  contour  irregularities  on  cavitation 
inception. 

The  committee  had  agreed  that  the  present  study  should  be  carried  out 
for  a  propeller  such  as  that  of  the  "Sydney  Express"  for  which  observations 
of  cavitation  have  been  made  on  the  full  scale  propeller  at  sea  and  at 
model  scale  in  various  laboratories.  However  it  was  found  that  the  pro¬ 
peller  fitted  to  the  "Sydney  Express"  has  blades  with  sharp  leading  edges. 
This  geometry  can  produce  laminar  separation  or  laminar  bubbles  on  the 
blade  leading  edges.  Since  anal/sis  methods  are  not  generally  available 
for  such  noncavitating  flows  it  was  decided  to  use  a  different  example  for 
the  present  calculations.  Accordingly,  the  skewed  propeller  fitted  to  one 
of  the  San  Clemente  OBO  class  merchant  ships  was  selected. 

The  design  parameters  for  this  propeller  are  documented  in  DTNSRDC 
report  P-500-H-04  [1,  2]*.  Only  a  few  of  the  characteristics  given  in  this 

* 

Numbers  in  brackets  indicate  references  cited  below. 
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source  are  required  for  the  present  calculations.  Moreover,  Reference 
[1]  summarizes  towing  tank  results  on  the  inflow  velocities  in  the  plane 
of  the  propeller  at  the  stem  of  the  hull.  Therefore  estimates  for  the 
relative  inflow  velocity  and  blade  element  attack  angles  are  also  available. 
These  data  provide  the  basic  ingredients  needed  for  the  present  work  in 
which  a  blade  element  at  the  .7  radius  point  is  taken  to  be  representative 
of  the  hydrodynamics  of  the  propeller.  Only  blade  surface  cavitation  is 
considered. 

PRELIMINARY  CALCULATIONS 

The  full  scale  0B9  propeller  has  a  radius  of  13  ft.  The  blade  chord 
at  .7  radius  is  6.557  ft.  The  towing  tank  model  has  a  radius  of  4.333  in. 
and  at  .7  radius  the  blade  chord  is  2.186  in.  The  factor  between  full 
scale  and  model  is  36. 

Reference  11]  gives  the  average  relative  inflow  angle  S.  as  11.72° 

at  the  .7  radius  station.  If  we  assume  a  ship  speed  of  17  knots  (28.39  fps) 

full  scale  and  3  knots  (=  5  fps)  for  towing  tank  tests  we  can  use  this 

information  to  estimate  the  Reynolds  number  range  for  the  present  calculations. 

If  we  base  the  Reynolds  number  on  blade  element  chord  length  and  relative 

•  5 

inflow  speed  we  find  for  the  model  that  Re  =  3  s  10  .  For  the  full  scale 

.  7 

propeller  we  find  Re  =  9  x  10  .  For  the  full  scale  ship  the  speed  of  17 
knots  is  close  to  its  maximum  speed.  Therefore  we  should  consider  a  range 
of  Reynolds  numbers  in  order  to  allow  for  lesser  speeds.  In  the  case  of 
the  model  we  need  to  consider  a  range  of  Reynolds  numbers  which  will  permit 
one  to  consider  towing  tank  tests  and  at  higher  speeds  cavitation  tunnel  tests. 
As  a  result  of  these  considerations  we  have  taken  the  following  range  of 


Reynolds  numbers  for  the  present  work. 
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Re  for  codel  tests  Re  for  full  scale 


3  s 

105 

io7 

io6 

3  x  10 

3  x 

106 

10s 

Turning  now  to  blade  element  geometry  at  the  -7  radius  point,  we  note 
that  a  NACA  series  6  casher  line  with  a  -0.8  was  combined  with  a  codified 
thickness  distribution  [3]  in  order  to  obtain  a  typical  blade  section.  The 
codified  profile  thickness  distribution  is  given  in  teres  of  percent  chord 
as  follows. 


x/c 

v/c 

x/c 

y/c 

x/c 

y/c 

0 

0 

20 

3.9957 

65 

4.3992 

.5 

-664S 

25 

4.3629 

70 

4.0545 

.75 

-S120 

30 

4.6382 

75 

3.6303 

1.25 

1.0436 

35 

4.8325 

SO 

3.1265 

2.5 

1.4663 

40 

4.9527 

85 

2.5437 

5. 

2.0655 

45 

5.000 

90 

1.8328 

7.5 

2.5253 

50 

4.9635 

95 

1.1453 

10. 

2.9073 

55 

4.8525 

97.5 

0.7485 

15. 

3.5213 

60 

4.6649 

100. 

0.3333 

According  to  Abbott  and 

von  Doenhoff 

[4],  the 

equation  of 

the  cean- 

line  for  series  6  profiles  is 


where 
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i  r  2.i .  is  .  li 

8  ■  -  PJ  [a  (2  ln  a  -  4>  +  4j  > 


"  '  £  I'1" 


a)2  ln  (1-a)  -  -|(l-a)2 


+  8 


The  ideal  attack  angle  is 


ot.  = 


Co  h 


i  2TT(a+l) 

where  h  is  the  maximum  camber  of  the  blade  section  as  a  fraction  of  its 
chord.  In  the  present  calculations  we  put  a  =  .8  and  eventually  combined 
the  mean  line  with  the  above  thickness  distribution  in  order  to  obtain  the 
profile  ordinates.  However  before  this  could  be  done  it  was  necessary  to 
convert  the  .7  radius  profile  for  the  0B0  propeller  into  an  equivalent 
two-dimensional  section  h-  removing  the  induction  effects  which  were 
necessary  because  of  lifting  surface  considerations  in  the  propeller 
design. 

On  page  40/,  Abbott  and  von  Doenhoff,  the  data  listed  for  the  mean 

line  having  a  =  0.8  show  that  if  C.  =  1.0,  a.  =  1.54°  and  the  camber  is 

*i  1 


=  .0679 


2-D 


From  Table  10  on  page  38  of  Reference  [1]  we  find  at  .7  radius  that 


f^|due  to  loading 


M 


=  k  =  2.025 
c 


2-D 
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at  Table  7  of  the  same  report  shows  that  — |  due  to  loading  is  .0323  at 
.7  radius.  Therefore, 


C  *  2D  2*°25 


=  .01595 


moreover,  — —  =  .0679  C0  from  which  it  follows  that 

2-D  i 


q  -  .2349 

A  .0679 


a±  «  1.54  x  .2349  =  .362' 


Continuing  with  data  from  Reference  [1]  we  find  from  Table  16  that  at  the 
.7  radius  station  the  variations  of  inflow  angle  from  towing  tank  tests 


+AB  “  5.08° 


-A8  =  -5.69' 


These  values  of  A(s  provide  for  the  distorted  inflow  caused  by  the  ship's 

wake.  These  values  can  be  summed  with  2  in  order  to  obtain  the  maximum 

and  minimum  values  of  attack  angle  for  the  two  dimensional  profile.  As  a 

result,  in  the  following  calculations  we  consider  three  attack  angles: 

a|  =  6.05°,  a.  =  0.362°,  and  a|  .  =  -4.72°. 

'max  l  'min 
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Having  determined  Cg  as  noted  above,  we  can  then  determine 

i 

the  profile  shape  by  adding  the  thickness  distribution  to  the  camber  line. 

A  plot  of  the  resulting  profile  is  shown  in  Figure  1.  Once  the  profile 
ordinates  have  been  defined,  one  can  find  the  pressure  distribution 
around  the  profile  in  accordance  with  potential  theory.  Three  pressure 
distributions  were  obtained  using  a  Douglas-Neumann  program  for  the  three 
values  of  a  determined  above.  The  resulting  pressure  distributions  are 
plotted  in  Figure  2  below. 

Once  the  pressure  distributions  are  determined  one  can  carry  out 

boundary  layer  calculations.  In  the  course  of  this  study  we  considered 

a  number  of  possibilities  for  executing  such  calculations.  We  have  chosen 

Truckenbrodt's  revised  method  [5]  because  of  its  simplicity.  Basic  steps 

in  this  calculation  and  the  final  results  are  given  in  Appendix  A  of  this 

report.  We  need  only  mention  here  that  such  calculations  were  carried 

out  for  upper  and  lower  surfaces  of  the  profile  of  Figure  1  for  the  three 

attack  angles  called  out  in  Figure  2  and  for  all  six  of  the  Reynolds 

numbers  listed  above.  In  all,  thirty-six  separate  calculations  were 

carried  out.  For  those  Reynolds  numbers  appropriate  to  model  tests  both 

laminar  and  turbulent  boundary  layers  are  expected  to  be  present  and  a  very 
£ 

simple  criterion  was  employed  to  fix  the  transition  point.  For  the  higher 
Reynolds  numbers  of  the  full  scale  propeller  the  boundary  layer  is  pre¬ 
dominantly  turbulent.  It  should  be  noted  that  high  precision  is  not  the 
aim  of  these  calculations.  Our  goal  is  to  provide  a  straightforward 
engineering  analysis  which  might  lead  to  useful  quality  control  requirements. 

* 

Appendix  A 
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EFFECTS  OF  ROUGHNESS 


Having  dispensed  with  the  necessary  preliminaries,  we  arc  free  to 
consider  the  matter  of  limits  on  roughness.  The  basic  reference  for  all 
that  follows  is  contained  in  the  J.S.R.  (1979)  paper  by  Arndt,  Holl,  Bohn, 
and  Bechtel  [6]  Figure  12  of  this  paper  shows  that  an  isolated  two- 
dimens'onal  triangular  roughness  element  will  produce  a  much  stronger 
effect  on  cavitation  onset  than  do  isolated  cones  or  triangular  irregularities. 
Moreover,  from  the  data  of  Figure  7  of  the  paper  it  appears  that  the  triangular 
2-D  element  produces  a  stronger  effect  than  the  circular-arc  2-D  element. 
Therefore  we  are  justified  in  considering  the  triangular  element  as  the 
dominant  case  to  consider  because  criteria  based  upon  it  will  certainly 
contain  other  forms  of  roughness.  We  shall  investigate  the  maximum 
roughness  height  which  will  permit  no  more  than  a  relative  effect,  E,  in 
increasing  the  onset  of  cavitation.  In  order  to  do  this  we  can  proceed 
as  follows. 

In  this  study  we  will  neglect  all  cavitation  scale  effects  and  put 
the  critical  cavitation  number  for  cavitation  onset  on  a  smooth  body 


equal  to  the  magnitude  of  C 


Thus 


mm 


a  =  -C 
cs  p 


min 


From  the  plot  of  Figure  2,  page  25,  below  we  see  that  at  the  three  values 
of  a  we  have 


a 


6.05° 

.362° 

-4.72° 


cs 

2.77 

0.35 

2.83 


x/ c 

0.01 

0.45 

0.01 
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The  critical  case  is  the  third  of  these  because  the  boundary  layer  will 
be  thin  and  the  critical  cavitation  number  is  large.  However  we  will 
consider  all  three  cases  in  order  to  compare  the  outcome  at  these  three 
attack  angles. 

Next  we  can  turn  to  the  results  summarized  in  Figure  7  of  Arndt  et  al. 
[6],  For  the  triangular  2-D  element  we  find 


where  is  the  onset  cavitation  number  for  a  protuberance  on  a  flat 
plate  having  no  pressure  gradient,  h  is  the  height  of  the  two-dimensional 
element,  6  is  the  boundary  layer  thickness,  U  is  the  velocity  at  the  edge 
of  the  boundary  layer  and  v  is  the  kinematic  viscosity.  Suppose  the  flat 
plate  has  the  same  length  as  the  profile  chord  c.  Then  the  preceding 
equation  can  be  expressed  as 


iq,  r,  .361/7*,  .165 
°fp  -  -152CR e)'196  $ 


Next  we  can  introduce  the  superposition  equation. 


a  D  =  -  C  +  (1  -  C  )  o, 
cR  Ps  ps  fp 


from  Equation  (6)  of  Arndt  et.al  [6].  Now  let  us  define  the  relative 
effect  of  roughness  by 


E  = 


o 

cs 


* 
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where  o  is  the  onset  cavitation  number  on  the  roughened  body,  Ocs  is 

the  same  quantity  on  the  smooth  body,  and  C  is  the  pressure  coefficient 

ps 

on  the  smooth  body.  Recalling  our  definition  for  E  given  above  we  can 
write 


E  =  .152 


(1-C  )  (Re)' 196 (h/c)' 361 

Ps 


(-C  )  (<5/c) 

Po 


.165 


Solving  for  h/c  we  get 


h 

c 


E2,77(6/C) 
.0054 (Re)' 


.457 

543~ 


2.77 


(A) 


In  the  calculations  to  follow  we  will  consider  two  values  of  E:  5 

2  77 

and  10%.  For  these  two  values  we  have  E  ’  /.0054  =  .046  and  .313 

respectively.  The  ratio  of  permissible  h/c  values  in  these  cases  is  6.8. 
The  following  table  gives  calculated  values  of  h/c  for  these  two  values 
of  E. 


-Re- 


a° 

3xl05 

106 

3xl06 

107 

3xl07 

00 

o 

6.05 

0.01 

4.36xl0"6 

1. 72xl0-6 

7.42x10" 7 

2.91xlO-7 

1 . 25xl0~7 

4.94x10  8 

0.362 

0.45 

1.69x10" 5 

9. 84x10" 6 

6.53xl0"6 

3. 22xlO~6 

1.61xl0"6 

-4.72 

0.01 

7.01x10" 6 

5.42xl0~6 

5. 69x10" 7 

1.75xl0”7 

7.48xl0~8 

3 .0xl0~8 

6.05 

0.01 

3. 19x10" 5 

1. 17xl0"5 

5.04xl0~6 

1.98xl0-6 

8.48xl0"7 

3. 36xlO-7 

0.362 

0.45 

1.15xl0"4 

6. 69xl0~5 

4.44xl0"5 

2.19xl0"5 

1.09xl0"5 

5.51xl0"6 

-4.72 

[ 

0.01 

4.81xl0"5 

3.68xl0“3 

3. 87xl0_6 

.  J 

1.18xl0"6 

5.08xl0-7 

2.04xl0-7 

model 

full  scale 
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The  values  tabulated  above  must  be  supplemented  by  additional  calculations 

for  the  minimum  value  of  C  on  the  lower  surface  at  a  =  0.362,  not  because 

P 

it  is  the  minimum  C  at  this  value  of  a  but  because  the  value  of  x/c  is  3.0% 
P 

and  although  =  -.159  at  this  point  the  boundary  layer  will  be  very  thin 
at  this  position. 


-Re- 

E 

3xl05 

io6  1 

3xl06 

io7 

mm 

io8 

5% 

1. 33x10”** 

6.97xl0”6 

2.99xl0-6 

1.18xl0”6 

l.llxlO”6 

5. 37xl0~7 

10% 

! 

9.07xl0~4 

4.74x10”*’ 

2. 03x10” 5 

8.02xl0"6 

7.55xl0”6 

3.65xl0“6 

model 

full  scale 

These  new  values  can  be  compared  with  appropriate  entries  in  the  table 
above  and  we  find  that  they  are  all  larger  than  those  associated  with  the 
minimum  pressure  coefficient.  Therefore,  it  appears  for  the  present 
three  cases  at  least  that  all  critical  roughnesses  are  those  calculated 
first. 

On  the  otherhand  it  must  be  recognized  that  the  data  which  underlie 
Equation  (A)  are  obtained  from  roughness  elements  which  were  not  smaller 
than  40  microns  ard  for  values  of  h/6  no  smaller  than  .0136.  This  value 
of  h/S  can  be  compared  with  the  displacement  and  momentum  thicknesses.  For 
a  laminar  boundary  layer  on  a  flat  plate  we  have 


^1  ^2 
~  =  .344  and  «  .133 

Data  in  Table  I  of  Reference  [ 6 J  for  triangular  elements  show  height 


ratios  as  small  as 
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^  =  .0136 


Since  the  experiments  were  carried  out  with  turbulent  boundary  layers, 
it  is  of  interest  to  compare  this  smallest  height  with  the  viscous 
sublayer  height,  6^.  We  can  estimate  this  quantity  by  using  (see 


Schlichtig,  op  cit.) 


r  =  5V 
*  '  V*  ’ 


where  v.  =  /r  /p.  Assuming  a  1/7  power  law  profile  and  allowing  for 
variations  in  pressure  coefficient  along  a  body,  we  would  have  for  the 
shearing  stress  at  the  wall 


-1/5 


p  '  V(1-Cp)C-0296)(Rec) 


After  some  manipulation  we  get 


‘  V/10/i^T 


If  we  put  0^  =  0  for  a  flat  plate  and  change  29  to  100,  neglect  (x/c)^^ 
and  replace  the  Re^^  with  Re  =  Uqc/v  we  get  the  admissible  roughness 
formula  in  Schlichting  (page  660-661).  Tt  appears  therefore  that  admissible 
roughness  can  be  as  large  as  four  times  the  viscous  sublayer  height. 

For  the  1/7  power  law  profile  we  have 
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6  _  -7  (x/c)4/5 

‘  Re1'5 


and  it  follows  that 


(x/c  Re)' 


where  we  have  retained  the  flat  plate  values.  For  a  Reynolds  number  of 
10^  and  x/c  =  1  we  get 


r  =  *005 


This  value  is  probably  the  largest  of  the  viscous  sublayer  heights 
encountered  by  Arndt  et  al.,  but  even  so  it  is  significantly  smaller  that 
their  smallest  h/6  noted  above.  If  we  were  to  multiply  the  ratio  6^/6 
by  4  we  find  that  h/6  is  in  the  range  of  admissible  roughness,  as  can  be 
seen  from  data  in  Table  12.3  of  Schlichting.  Since  the  present  analysis 
appears  to  show  a  possible  effect  on  cavitation  of  roughness  elements 
which  are  smaller  than  those  of  the  experiments  it  would  be  of  great 
interest  to  extend  these  experiments  to  even  smaller  roughness  sizes.  When 
one  couples  this  finding  with  the  fact  in  many  situations  encountered 
in  the  tawing  tank  or  water  tunnel  that  laminar  boundary  layers  are  present, 
the  extension  of  flat  plate  roughness  data  to  include  cavitation  in  laminar 
boundary  1  avers  would  also  be  beneficial. 
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At  this  time  it  appears  that  our  attempt,  using  2-D  triangular  roughness, 
to  define  criteria  for  the  onset  of  cavitation  on  marine  propellers  is  hampered 
by  a  lack  of  data  although  the  methods  by  which  such  rational  standards  might 
be  developed  seem  straightforward. 

As  we  have  noted,  the  preceding  calculations  assume  an  isolated  two- 
dimensional  roughness  of  triangular  shape  because  of  its  strong  effect  on 
cavitation  inception.  On  the  other  hand  if  we  x*ere  to  go  to  the  other 
extreme  and  consider  distributed  roughness  only,  we  may  be  justified  in 
using  the  usual  criterion  based  on  skin  friction  which  defines  admissible 
roughness.  We  base  this  assertion  upon  the  correlation  given  in  Figure  6 
of  Reference  [6].  In  this  case  we  would  have 


h  100 

Re/l-C  (x/c  j 
P 


(B) 


for  the  admissible  height  of  distributed  roughness.  In  this 
would  be  no  measurable  effect  on  cavitation  inception;  that 
of  the  above  terminology  E  =  0.  Clearly,  the  critical  value 
be  found  at  C  I 


P 


min 


case  there 
is ,  in  terms 
of  h  would 


EFFECTS  OF  A  CONTOUR  IRREGULARITY 

The  formulation  of  a  method  by  which  one  can  examine  the  effects  of 
variations  in  profile  contour  seem  somewhat  less  demanding  than  the 
preceding  considerations.  One  way  is  to  use  potential  theory  to  calculate 
the  pressure  distribution  of  the  deformed  profile  from  which  the  cavitation 
performance  can  be  reassessed.  The  chief  problem  with  this  approach  is 
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that  not  all  organizations  have  a  computer  program  possessing  the 
"numerical  flexibility"  which  would  enable  them  to  handle  small  bumps 
or  waves  without  difficulty.  For  those  who  do  not  have  such  a  "flexible" 
program,  we  would  propose  the  following  approximate  procedure. 

The  basic  idea  here  is  to  use  the  superposition  equation  from 
Reference  [6]  as  quoted  above.  For  O ^  one  would  simply  use  the  magnitude 


of  C 


of  one  of  the  bumps  given  in  Appendix  B  of  this  report.  See 


min 


Figure  3  below.  The  conditions  under  which  this  procedure  is  valid  are 
shown  on  Figure  2  of  Reference  [6].  The  curves  and  the  data  given  there 
show  that  if 


h/S  >^4  , 

then  the  proposed  procedure  should  work  satisfactorily.  Here  h  is  the 
maximum  height  of  the  bump.  Otherwise  the  methods  which  apply  to  smaller 
roughness  elements  as  discussed  above  should  be  used. 

In  terms  of  the  bumps  and  pressure  distributions  of  Appendix  B  we 
may  note  that  the  parameter  e  is  given  by 

e  =  h/£  , 

where  £.  is  the  length  of  the  bump  (taken  as  unity  in  those  calculations) 
and  h  is  the  actual  height  of  the  bump.  Therefore  we  can  write  the 
criterion  of  applicability  as 


In  those  cases  where  the  observed  bump  is  not  well  represented  by  either  of 
those  given  in  Appendix  B  one  can  develop  the  measured  contour  deviation  about 
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a  straight  line.  Then  he  can  superpose  bumps  having  various  values  of  £  along  this 

line  until  the  observed  deviation  from  the  line  is  fairly  well  represented.  After 

that  he  can  superpose  the  pressure  distributions  from  all  bumps.  Values  of 

x  (or  x/Jt),  n/e  and  C^/e  tabulated  in  Appendix  E  could  be  useful  in  this  phase 

of  the  calculations.  The  0^  will  be  equal  to  the  magnitude  of  C  j  obtained 

Hmin 

from  the  composite  pressure  distribution. 

As  a  simple  illustration,  let  us  suppose  that  only  one  bump  is  needed  to  . 
characterize  at  least  approximately  a  typical  deformation.  If  we  settle  on  Case  ?2 
from  Appendix  B  we  see  that 


Then  since 


f 


=  -.68e 


mm 


we  have 


E  =  .68 


e 


Now  if  we  put  E  =  . 1  we  have 


e 


.  147 


Finally  we  can  use  the  "criterion  of  applicability"  from  above  and  write 
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j  <  -0368 


_r 


1-C 


and  ve  also  see  that  it  leads  to  the  alternate  expression 


e  <  -147 


-C 


1-C 


As  an  application  of  these  results  suppose  that  we  consider  the 
critical  case  fron  the  plot  of  Figure  2  which  is  found  at  a  =  -4.72- 

C  I  for  this  case  is  -2-84  and  we  find 

P  - 
jnin 

e  <  .109 


f<_  .0272 


On  the  ether  hand  we  find  fron  Appendix  A  that 


~  =  .00013  when  Re  =  10** 


Hence  it  follows  that 


.005 


These  data,  together  with  the  fact  that  the  buap  is  given  by 


2 

rj(x)  =  £(1-4-0 
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where  x  is  measured  from  the  center  of  the  hunp  and  is  an  actual  distance 
normalized  by  f.  and  n  is  the  normalized  ordinate  of  the  bump;  serves  to 
characterize  a  permissible  deformation  which  will  produce  a  relative  error 
in  cavitation  number  of  about  102-  The  effect  of  more  complex  bumps  can  be 
explored  by  use  of  the  composite  pressure  distribution  discussed  above. 

CONCLUSIONS 

An  investigation  which  makes  use  of  the  most  complete  information 
available  on  the  effects  of  surface  asperities  with  respect  to  cavitation 
inception,  has  sought  to  determine  how  these  data  might  be  used  to  formulate 
quality  control  criteria  for  cavitation  inception  on  model  propeller  blades. 
The  OBO  propeller  was  used  as  an  example.  It  was  found  that  the  dominant 
effect  is  caused  by  the  two-dimensional  isolated  triangular  element.  In 
this  case  the  data  are  not  complete  enough  to  permit  the  formulation  of 
a  roughness  criterion.  On  the  other  hand,  if  one  considers  the  effaces  of 
distributed  roughness  or  single  bumps  which  protrude  well  outside  the 
boundary  layer  one  can  formulate  methods  for  obtaining  simple  quality 
control  criteria  for  these  two  cases. 
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PRESSURE  COEFFICIENT 


CURVES  SYMMETRICAL  ABOUT  VERTICAL 


Figure  3.  Pressure  Distributions  from  Appendix  B  of  Two  Bumps  on  a  Flat  Plate. 

Case  //I :  Linearized  "Circular"  Bump  (Parabolic  Contour).  Case  112: 

•  Bump  with  Cusp. 
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APPENDTCIES 


A:  Boundary  Laver  Thickness 

Truckenbrodt' s  Revised  Integral  Method. 

See:  Boundary-Laver  Theory  by  H.  Schlichting,  7th  Edition,  McGraw-Hill 
Book  Company,  New  York,  1980,  Chapter  XXII. 

Basic  Equations  Used:  Let  x  be  arc  length.  Then 


<5~U  (x) 

&3  (x)  =  — ~ 


c 


c  U 


ReA-C  (x) 

P 


where 


6 

6^(x)  =  energy  thickness  =  j  [1- (u/U (x))^]  (u/U)dy  , 

o 


V  =  kinematic  viscosity, 
c  =  profile  chord  length, 

U(x)  =  velocity  at  edge  of  boundary  layer, 

u(y)  =  velocity  profile  at  any  arc  length,  x, 

U  =  free  stream  velocity 
o 

Re  =  Reynolds  number  based  on  free  stream  velocity  and  chord  length 
6(x)  =  boundary  layer  thickness  at  the  arc  length  x 
Laminar  Boundary  Layer: 

If  we  assume  that  u(y)  can  be  approximated  by  the  Blasius  profile  we  get 


-  =  4.788  — 
c  c  ’ 


as  can  be  found  from  a  simple  Simpson’s  rule  integration  as  shown 


below. 
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Calculations  for  the  ratio  6/6  for  the  Laminar  Boundary  Layer. 

s 

Blasius  Flat  Plate  Solution  (see  Schlichting  p.  139,  7th  Edition). 


63  ‘  5 


no 

f*  (l-[f'  ]2)d£  =  1 1 


C 

F 

K 

F 

1. 

0 

0 

22. 

4.2 

.06284 

2. 

.2 

.06612 

23. 

4.4 

.04653 

3. 

.4 

.13043 

24. 

4.6 

.03373 

4. 

.6 

.19107 

25. 

4.8 

.02397 

5. 

.8 

.24616 

26. 

5.0 

.01669 

6. 

1.0 

.29392 

27. 

5.2 

.01140 

7. 

1.2 

.33272 

28. 

5.4 

.00764 

8. 

1.4 

.36128 

29. 

5.6 

.00502 

9. 

1.6 

.37876 

30. 

5.8 

.00323 

10. 

1.8 

.38489 

31. 

6.0 

.00204 

11. 

2.0 

.38000 

32. 

6.2 

.00126 

12. 

2.2 

.36505 

33. 

6.4 

.00078 

13. 

2.4 

.34158 

34. 

6.6 

.00046 

14. 

2.6 

.31154 

35. 

6.8 

.00026 

15. 

2.8 

.27708 

36. 

7.0 

.00016 

16. 

3.0 

.24045 

37. 

7.2 

.00008 

17. 

3.2 

.20366 

38. 

7.4 

.00004 

18. 

3.4 

.16846 

39. 

7.6 

.00002 

19. 

3.6 

.13616 

40. 

7.8 

0 

20. 

3.8 

.10756 

41. 

8.0 

0 

21. 

4.0 

.08311 

->■ 


im 


•  % 


n  =  40 
h  =  0.2 

Using  a  TI  59  with  IfL  10  for  Simpson’s  P.ule  we  get  I  =  1.04434.  Therefore 


56. 


6  = 


1^34  •  '"7877153 


The  calculator  results  are  shown  on  the  next  page. 


n  \Cn 
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4.788 

/Re 


ft 

[l-c  (C)l  dC 

o 

1/2 

_  _5_ 

|\l-C  (£)l2'5d€ 

J  P 

o 

.917  u-ye,)]3 

✓Re 

U-Cp^)]3 

(A-l) 


It  is  worth  observing  that  when  C  =  0  in  the  laminar  boundary  Equation 

P 

(A-l)  that  we  get 


—  =  — f— ] 

c  /Ri  W 


1/2 


which  is  the  well  known  Blasius  result.  Moreover  in  Equation  (A-l)  the 
quantity  £  is  given  by 

K  =  x/c  < 

and  as  is  customary,  at  the  forward  stagnation  point  one  puts  6=0. 

Equation  (A-l)  is  used  to  calculate  the  laminar  boundary  layer  thickness 
from  the  nose  of  the  profile  to  the  transition  point.  For  many  purposes  it 
is  permissible  to  use  a  simple  cricerion  in  u.Jer  to  estimate  the  location 
of  transition.  One  can  suppose  that  transition  may  occur  at  an  arc-length 
Reynolds  number  of 


Rec  =  3  x  10' 


or  in  the  neighborhood  of  C  ,  whichever  occurs  first.  The  criterion 

^min 

is  probably  most  appropriate  if  the  outer  flow  is  highly  turbulent  as  in 
the  wake  of  a  ship  although  its  chief  virtue  is  simolicity  rather  than 
accuracy. 

Once  the  transition  point  :<1  has  been  selected  one  can  calculate 
the  laminar  energy-layer  Reynolds  number. 
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S(x.)  Re/l-C  (x.) 

R  (x  )  =  _ i _ P  -1 

V  r  c  4.788 


(A- 2) 


for  the  transition  point  only.  One  assumes  that  the  energy  layer  6^  is 
continuous  at  transition. 


Turbulent  Boundary  Layer 

Once  the  laminar  boundary  layer  thickness  distribution  and  the  location 
of  transition  have  been  determined,  calculations  for  the  turbulent  boundary 
layer  can  be  made.  The  formula  for  the  energy  thickness  Reynolds  number  is 


j  [i-cp(C)]1,65dC 


1 

1.152 


„  ,  „  r„  ,  ,  ,1.152  .  Re  S1 

Ro(x)  =  [R-(x.)]  +— - 

J  J  1  v  [1-C  (x)T*x:> 


(A- 3) 


where  R^(x^)  is  given  by  Equation  (A-2),  E,  =  x/c  as  before  and  £  =  x^/c. 

Schlichting  suggests  that  for  v'  one  should  use  v'  =  80v.  The  boundary 
layer  thickness  can  be  determined  from  the  energy  thickness  if  it  is 


assumed  that  the  velocity  profile  can  be  approximated  by  a  one-seventh 
power  law  profile  which  is  known  to  give  rather  good  results  for  flat- 
plate  skin  friction  at  Reynolds  numbers  less  than  10^.  As  a  result  one 
finds  that 


6  =  40  63(x)/7  , 


from  which  it  follows  that 
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6(x)  _  5"71A  R3(x) 


ReA-C  (x) 
P 


(A-4  ) 


In  the  course  of  preliminary  calculations  based  upon  Equations  (A-l) 
through  (A-4)  the  recommended  value  of  v*  =  80  v  was  used.  Surprisingly, 
when  laminar  and  turbulent  boundary  layer  calculations  were  started  from 
the  nose  of  a  profile  it  was  found  that  the  laminar  boundary  layer  thick¬ 
ness  was  larger  than  the  turbulent  boundary  layer  thickness  for  a 
considerable  range  of  arc  lengths  starting  from  the  nose. 

Therefore  we  conducted  a  short  investigation  to  see  if  a  better  choice 
for  the  value  of  v'  could  be  made.  We  considered  a  flat  plate  of  length  c 
at  Re  =  10^.  The  laminar  boundary  layer  thickness  for  this  case  can  be 
found  from  the  formula  at  the  top  of  page  30-  The  corresponding  turbulent 
boundary  layer  thickness  for  a  one-seventh  power  low  profile  is  known  to  be 


6(x)  _  .37 (x/c) 

c  D  1/5 

Re 


4/5 


and  if  one  puts  C^(£)  =  R^(x^)  =  0  in  Equations  (A-3)  and  (A-4)  it  follows 
that 


600 


5.714  f— -1 

lv'  C> 


.868 


(Re) 


.132 


If  one  equates  these  two  values  of  turbulent  boundary  layer  thickness  he 
finds  that  the  ratio  v/v'  is 
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V 


.0427 


v'  [(x/c)Re] 


.0452 


(A-5) 


Equation  (A-5)  can  be  used  to  explore  various  possible  values  for  v'.  For 
example  we  have  used  Re  =  10^  as  noted  above  and  we  considered  a  few 
values  of  x/c  as  tabulated  below. 


Re  =  10 


>c 

c 


V  1 

v 


.5  65 
.1  58 

.05  38 


Table  showing  the  effect  upon  possible 
values  of  V’  for  three  transition  point 
locations  on  a  flat  plate. 


Next,  two  likely  values  of  v'  were  used  to  compare  1/7  power  and  Trucken- 
brodt  boundary  layer  thicknesses  directly  at  Re  =  10^  and  x/c  =  .08  in 
which  case  6/cl^^  =  .00309. 


For  V  =  50v  ,  — 


=  .00345 


For  v'  =  55v 


=  .0032 


As  a  result  of  these  comparisons  it  seemed  best  to  take 


V'  =  oov 


for  the  present  calculations.  Nevertheless  as  a  final  illustration  of 
the  effect  of  changing  from  v'  =  55v  to  v*  =  80v  on  a  typical  profile  is 

C 

illustrated  by  the  data  on  page  39  belo’w,  which  are  given  for  Re  =  3  x  10' 
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and  Re  -  10  when  v*  =  80v.  These  data  can  be  compared  with  the  corresponding 
results  for  v'  =  55v  on  pages  36  and  38  respectively. 

Numerical  Results  and  Calculator  Program 

Equations  (A-l)  through  (A-4)  were  coded  for  execution  on  a  Tl-59 
programmable  calculator.  For  the  record,  we  give  a  description  of  the 
procedure  and  other  details  of  this  phase  of  the  calculations. 


Data  Registers 


Roo  - 

for  DSZ  index 

R61  - 

blank 

R01  - 

used ,  R^j 

and  R^  laminar 

R62  - 

[1-C  ],1  <  j  <  51 

Pj 
j 

R03  - 

h  in  Simpson's  rule 

R63  - 

R04  - 

I,  Simpson' 

's  rule  result 

R64 

k 

i 

in 

O 

j  .  =  total  no  of  C  entries 

IH3X  p 

R65  - 

Re 

R06  - 

f0’  k  =  6 

R66  " 

II  from  R„, 

R07  ' 

fr  k  =  7 

•  in  Simpson's  rule 

R67  " 

6/c  laminar 

R08 

f2,  k  =  8  . 

R68  - 

R^(x)  turbulent 

R09  - 

jC^)  =  an 

odd  integer  denoting 

R69  - 

f^.,  6  <  k  <  8  where 

R10  = 


location  of  (transition) 

C  (0)  -  R.  =  C  (x  ),  a  total 
p  60  p  max 


fk  -  1>V 


of  51  C  (x 
P 

1  ''  j  <  51 


C  (x)  entries  can  be  used 
P 


Preliminaries 

(a)  Partition  memory  j7|  ]2^d|  |o£|  |l7j,  display  399.69.  To  check, 

i2^di  m  mi. 

(b)  Enter  program  from  magnetic  card  Joj  ilNVj  find)  jWF.ITE j  for  both  1  and  2. 


. . . . . . . 
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Procedure 

(a)  entar  data:  ®  ISl  ISI  1^1  IS®  !°SI  l3SII  iMI  I® 

mi  suit  mi  _ 

(b)  entar  CPj:  ISl  153  ISl  ©  ISl  ®  ^ 

J  R/S  |  this  moves  pointer  to  R^q*  I^PjJ  l^£_l 

lgg  __ 

(c)  Run  program:  |CLR|  |RST|  |R/s| 

- H - 

Program  Listing  and  Data  Register  contents  given  on  pages  55  -  57.  Note 
value  of  V’  in  Program  Memory  addresses  306  and  307. 

Specific  boundary  layer  calculation  results  for  GBO  propeller  are  given  on 
pages  36  ~  5A. 


^llillllPliiiiulliUllHklilllliiillil'ItillilillltlUl  allill  iliSlIli'iillhiJlhiililiDlllllllillllltllllliritlillllllllilt1'  illilllllllllllllllllllllllllllilllll 


L 


i  J ,  h 
■Vi  fi;‘ " 


0,  64 
, 0234323377 


0.  .:8 


Upp.VL 


tt'S.O'? 


Ml  i 1 


nmnpi:-:-?7l9 

49 


!l  Ih,  1  I  16 


t 
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Ur  hr 

.  0182840755 

t  U  i.  J  6  4675 8 8 2 

LL  T  b 

0.  44 

,  0196251926 

« 0092221973 

u  i  o 

Ur  -rO 

.  n  ’  i  f=i2i  4  f, 

»  U  1  ■  9  9  6' 1 1 2 0  4 6 

U  7  1 U  9 1 


i  uooonno! 


0.  36 
0030 1 42365 


«  UUUU 


O-"' 


0.  ?s 

=  0202417  9  £  3 


0.  24 
0  =+  i  ■  6 


j"_  j" j  ‘  J.  -  ' ' 


.  i  "»0 
■T  1  r**i 


Or  32 


Lcrurt^c. 


L*  I 


B  -^7”:^ ' : 9DS  m  S  \  E 1 1 
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B:  Perturbed  Pressure  Distributions  Caused  by  Surface  Irregularities 


In  this  discussion  we  include  two  bumps  because  they  lead  to  the 
simplest  computations.  The  bumps  presented  are: 


1st  case: 

n(x) 

=  e(l  - 

4  x2) 

,  I x I  <  1/2 

* 

O 

II 

« * 

| x |  >  1/2  . 

2rd  case: 

n(x) 

=  e(l  - 

4  x2)2 

,  |x|  <  1/2 

;  =o 

| x |  >  1/2  . 

For  Case 

1-. 

we 

have 

n'(x) 

= 

-8ex 

.  1 

x|  £  1/2 

;  n'(x)  =  0 

elsewhere 

* 

n"(x) 

= 

-8e 

.  1 

x|  <  1/2 

;  n”(x)  =  0 

elsewhere 

• 

For  Case 

2; 

:  we 

have 

n’(x) 

= 

-16ex(l  -  4 

x2) 

|x|  <  1/2  ; 

n’(x)  =  o 

elsewhere 

n"(x) 

= 

-16e 

(1  -  12 

x2)  , 

I x I  <  1/2  ; 

n"(x)  =  o 

elsewhere 

Note  that  both  bumps  have  discontinuous  curvature  at  Jx{  =  1/2  but  that 
in  Case  2  the  slope  is  continuous.  We  have  considered  bumps  having  a 
continuous  slope  and  curvature  at  |x|  -  1/2  but  the  series  convergence 
is  not  too  good. 

Case  1 

Consider  the  bump  for  which  n’ (x)  =  -8ex  as  noted  above. 


2  =«MH{  ♦  l/{  i 
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In  accordance  with  the  ideas  of  thin  airfoil  theory  we  consider  the  complex 
velocity,  W(z)  =  U(1  +  u  -  iv),  to  be  invariant  at  corresponding  points  in 
the  z  and  4  planes.  Next  we  consider  the  mapping, 

z  =  (x  +  iy)  =  |  (4  +  -|)  . 

i  0 

on  the  unit  circle,  4  =  e  ,  and  2x  =  cos9.  Accordingly,  n’(x)  =  -4ecos6. 
But  in  accordance  with  the  scheme  of  linearization  the  boundary  condition, 
v(x)/U  =  n'(x),  is  to  be  evaluated  at  y  =  0  on  the  real  axis.  Thus  we 
have 

n'(x)  ,  | x J  _<  on  the  bump  , 

0  ,  |x|  y  off  the  bump 

But  we  know  that  for  flows  without  circulation 

J(z)  =  U[1  +  w(z)]  can  be  written  with 


/  %  t*  3 .  lb . 

w(z)  =  L  _JL _ 1 

1=1  c1 


(where  w  -  u  -  iv  is  the  perturbation  velocity) 


in  order  that  w(z)  0  as  (  »  and  that  the  flow  direction  at  «  will  be 

invariant  under  the  mapping.  This  latter  aspect  is  accounted  for  in  the 
Joukowski  transformation. 

i6 

In  general  terms,  4  =  re  v.  But  on  the  unit  circle 


4  =  e 
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Therefore,  the  real  axis  in  the  interval,  (y  =  0,  | x |  <_  1/2)  is  on  the 
unit  circle  with  0  8  _<  tt.  The  point  (y  =  0,  x  =  0)  and  at  the  center 

of  the  bump,  where  n'(0)  =  0  and  n(0)  =  e,  maps  into  the  point  6  =  tt/2,  r  =  1. 
Accordingly,  we  have 

00  oo 

w  =  u  -  iv=  E  (a.  +  ib.)e  =  E  (a.  cos  i8  +  b.  sin  i6) 

..xx  .  .  x  1 

i=l  i=l 


OO 

-  i  E  (a.  sin  i0  -  b.  cos  i0) 
1  1 
i=l  X 


Applying  the  boundary  conditions  on  the  bump  we  have 


-iv  =  -i  E  (a.^  sin  i0  -  b^  cos  i6)  =  +i8ex,  |x|  _<  y  =  0,  |x|  >  — 
i=l  1  1 

Therefore  in  the  £  plane  we  have 


v 


E  (a.  sin  i8  -  b.  cos  i0)  =  -Ae  cos6  on  £  =  eX^ 
i=l  1  1 


But  if  v  *  0  at  points  on  the  real  axis  outside  the  unit  circle  and  on  |s|  =  1 
if  v(-9)  =  -v(9) ,  in  order  that  v(9)  is  odd  with  respect  to  8  so  that  it 
will  represent  a  thickness  distribution  and  not  a  camber  function,  we 
must  consider  the  expansion  of  cos6  as  a  sine  series.  In  this  case 
b.  =  0  and 

l 


(cos8)  sin  iG  d0 


-16si 

*(i2-l) 

’ 

0 


1  2y  A  *  6,  *  •  •  1 


otherwise 
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As  a  result  of  this  finding  it  turns  out  that  the  u  component  is 


given  by 


16e  „  i  cos  i0  ^ 

u  - - Z  - s -  ,  (1  is  even  onlv) 

V  i=2  (i  -1) 


and  since  the  linearized  bernoulli  equation  is  =  -2u,  we  have 


„  32e  t-  i  cos  i0  i0 

C  = -  Z  - ^ -  ;  on  C  =  e  ,  and  i  even 

p  71  i=2. . .  i2-l 


At  this  point  we  need  to  pause  and  consider  the  convergence  of  0^(8). 

In  accordance  with  the  Laurent  expansion  for  w  given  on  page  59  it  is  clear 

i0 

that  for  points  outside  and  on  the  circle  we  must  have  t,  =  re  and 


C  (£)  =  I  ^TC°y  ^  ,  with  i  being  even 

P  "  i=2  rx(i  -1) 


But,  on  the  circle  r  =  1  and  on  the  real  axis  r  1  and  9  =  0,  u.  In 
particular,  if  r  >  1  the  series  converges  for  all  0.  When  r  =  1  the 
series  diverges  at  0=0,  n.  This  can  be  verified  by  employing  the  usual 
tests  for  convergence  of  power  series.  It  can  also  be  seen  for 
sufficiently  large  values  of  i  when  0=0  (or  u)  that  the  series  becomes 
equivalent  to  the  harmonic  series  which  is  known  to  be  divergent.  There¬ 
fore  we  see  that  in  terms  of  thin  airfoil  theory  we  can  consider  the 
corner  between  the  wall  and  the  bump  as  a  linearized  stagnation  point. 

Since  the  series  at  other  points  is  convergent,  we  can  compute  values  of 
Cp  with  it,  but  we  should  not  expect  rapid  convergence.  Indeed,  as  we 
shall  see  below  even  though  we  take  many  terms  (as  high  as  40)  the  accuracy 
is  not  great. 
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Case  2:  n'(x)  =  -16  ex  (1  -  4x2)  =  -8  e  cos  8  sin2  8 

2 

-  -4  e  cos  8  sin  8  =  -2  e  (cos  8  -  cos  38) 

Note  that  this  case  contains  the  first  case  as  part  of  the  solution. 
Thus  we  see  that  the  term  -2e  cos0  will  lead  to  a  contribution  for  u(G) 
which  is 


GO 

-8£  z 

n 


1  cos 


i0 


i  =  2,  4, 


i=2  (i  -1) 


It  remains  for  us  to  evaluate  the  term  2e  cos  30.  In  this  instance  we 
have  to  consider 


cos  30  sin  iO  dG  =  ( 


8ei 


« (i2-9) 


n  =  2,  4,  6,  . . .  •  , 

otherwise 
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If  we  add  these  results,  terra  by  terra,  we  obtain 


8c 

i  i  ' 

i3-i  -  i3+9i 

1 

I  64c  i 

71 

o  2 

i  -9  i  -1 

n 

(i2-l)(i2-9)_ 

11  (i2-l)  (i2-9) 

Therefore  we  have 


„  .  .  128c  _ 

C  (x)  =  -  — —  I 
P  “ 


i  cos  i9 


i  2  2 

i=2  r(i  -1) (i  -9) 


32e  y  i  cos  ie  4 
11  i=2  r1(i2-l)  19-i2 


where  i  =  2,  4,  6  ...  ,  as  before. 

Note  that  we  have  written  this  expression  so  that  each  terra  of  the  sum 

for  0^  can  be  obtained  from  the  corresponding  terra  of  the  first  case  by 

2 

multiplying  it  by  the  factor  4/(9-i  ). 

The  convergence  of  this  series  is  no  problem.  As  i  “  each  term 
2 

approaches  1/i  ,  and  it  is  known  that  this  will  produce  a  uniformly  and 
absolutely  convergent  series  both  on  and  off  the  circle  r  =  1. 

Similarly,  it  is  worth  noting  that  the  formulae  for  n(x)  in  the  two 
cases  are  very  simply  related.  In  fact,  if  we  want  to  get  case  2  from 
case  1  vc  need  only  take  q (x)/e  for  case  1  and  square  it.  As  a  result  it  is 
fairly  easy  to  program  a  routine  which  will  do  both  cases  in  sequence  and 
this  has  been  done  below. 

It  should  be  noted  that  the  2nd  case  was  computed  on  the  basis  that 

2  2 

the  multiplicative  factor  4/(9-i  )  was  2/(9-i  ),  and  this  is  incorrect. 

Therefore  we  should  multiply  all  results  for  C  (x)/c  which  are  tabulated 

P 

below  by  a  factor  of  2.  This  error  was  corrected  in  the  case  N  =  40. 
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Notes  on  T.I.  -  59  program  for  finding  C^(x")/e  for  the  tvo  bumps: 


Case  1:  n(x)  =  e(l  -  4  x  )  , 

Case  2:  n(x)  =  e(l  -  4  x2)2 
Procedure:  Enter  program  from  keyboard.  See  below. 

Enter  Input  Data:  Data  Register  Address 

Total  number  of  terms  in  sum  =  N 
Increment  between  successive  x  =  Ax 


R01  *  |ST0!  01 


Maximum  value  of  x  =  x 


max 


r07  IMI  07 
ro8  -  !Mi  os 


In  order  to  execute  program  after  input  data  are  entered: 
Press:  jCLRj 

[rst| 

\m\ 


If  it  is  desired  to  run  Case  2  only,  enter  input  data  as  before.  Then 
Press:  |CLR[ 

liEIl  I'St  fig  I 
iM!  IM1 
IlZIi 


Storage  Map: 


R00  for  DSZ(n) 

s04  for  rW 

R0S  f°r  S | max 

\l  for  8 

R05  f"  1 

R09  f°r  Cp(K)£ 

R02  for  X 

R06  f°r  n^"^£ 

R10  for  (x-1/2) 

R03  f°r  6 

>r 

R„_  for  Ax 

0/ 

used  in 

summation  of 

Cp 

* 

Input  quantities. 

- H - 

Data  calculated  for  two  cases,  N  =  20  and  N  =  40  are  given  below.  These 
are  followed  by  the  program  listing  and  other  information.  A  plot  of 
Case  1  and  Case  2  is  given  in  Figure  3,  page  26  above. 


m 


& 

1 


ii 

•< 


m 


i 
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i.  4 

0. 

i.  £347862 


1,  45 
0. 

<  2176112402 

1,  5 

0. 

,  2022962247 

i,  55 
0. 

. 1885765801 

1  •  6 

0. 

.  1 1"  6 2 3 3 9944 

1.  65 
0. 

1  £  f| ■"« £  C* C1  *? n 
u  JL  ©  L  o  ©  O  O  1  O 

lx  7 
0. 

=  1549831282 

1=  75 
0. 

. 1457944004 

1.  8 

0. 

. 13741 18072 


i «  o  ■  j 

0. 

0 1  1297424 U 2 


0. 

0. 122706431 

1.  95 

0. 

0,  1 1  62351  22 


0. 

1 102686982 
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0.  7 

0. 

lr  294791052 

Or  75 
0. 

i  lr 054641806 

Or  8 

0. 

. 8814365151 

0.  85 
0. 

Or  750894088 


U. 

. 6492842871 

0.  95 

0. 

. 5682095794 


ir 

0. 

r  5022282621 
lr  05 

Or 

.  4476533426 

1=  1 
f  j 

0=4  u  1  y  “  r.  *4  ^  4 


u, 

36308634 , 4 


1 5  £5 


r«  i  -  n  *=:  t  i 

*  *-•  j.  JL  U  •  j  i  1 


i  a 
0 . 

.  2  6  n  0  5  5  0 :  5 


1 


imnjiN'WW' 


_ ' _ _ _ _ . - 


})ATfi  iZ£4tST£{ZS  | 

17.  /fen-bSi:  00 

40/  N  oi 

0. 35  X  02 

7953933302  0  03 

1.  a.  04 

3.  X  05 

0.  2601  06 


2601  r\(,%)/e  06 

0=  05  AX  07 


2.  Xyno^y.  03 

0664347757^ C A/ 6  09 
-0.  1 5  (vc-  '/iO  1 0 
- .  03  3 3 3 3 3 3 3 3  use’-?  1 1 

£co*is.*U >0*  ^iccf  J2 

/"  0.  ytO-"^C  1  - 

na-t  0.  14 

imccC  0.  15 

v\  0.  16 

-0.  17 

0.  IS 

0.  >*19 


0,.  ddcCrt^  20 

0.  21 


L  ctJtb<C^-  UjkcdL 
OP  ©2 


13  C 

15  E 

16  fl! 
43  EXC 

17  B* 
IS  C1 
19  D* 

14  D 

1 1  ft 

12  B 
10  E! 
66  PflU 


P/Ur^i Avn. 


000 

o 

l_  t— 

I NV 

4 

*  i 
#< 

no  l 
002 

S6 

01 

STF 

01 

i  f ' 
i 

B 

A 

" 

003 

00 

U 

p.vV' 

J 

004 

CL  > 

•i*70 

X 

V 

'  f 

005 

02 

02 

< 

-i 

*r 

s 

n  n  c 
>-•  <„•  >-* 

007 

003 

43 

01 

42 

Di-i 

I  S 

01 

ST0 

IfV 

pse 

M 

?V~s 

010 

00 

0 

/<-v  1 

on 

42 

ST0 

Or 

Cf  , 

012 

09 

QC| 

013 

02 

z 

irr  i 

014 

42 

STD 

•  A  ‘ 

* 

015 

05 

05 

016 

43 

RCL 

017 

02 

02 

013 

*“?  tr 
i  -J 

- 

>-y 

0 1 9 

93 

.  c> 

020 

05 

5 

kn- 

021 

95 

= 

v  .  • 

+**t 

022 

42 

STD 

023 
- - 

10 

10 

_ 

1  -■ 

14 

15 

16 
17 
IS 

^19 


«  U  O  1 

'L+  032 


4  ^cc034 

ytc  035 

036 

1-  037 

033 

_ 039 

MoIa-v.  U  4  0 

i\c<p  o  vV*  0  4 1 
042 
043 
044. 
4 

ruA 


047 

043 

049 

43 

06 

99 
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71 

IS 
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,~i  cr  v 
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£A  '-J 

09 

99 
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I^aau 

PRi^f; 

055 

93 

HDV 

056 

43 

RCL  , 

057 

02 

lj>* 

058 

85 

* 

059 

060 

4  ; 

^r-I* 

07 

n  !-1-  a 

07 
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